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^ Abstract 

A diagrammatic approach to calculate n-point correlators of the primordial curvature perturbation £ was 
developed a few years ago following the spirit of the Feynman rules in Quantum Field Theory. The method- 
ology is very useful and time-saving, as it is for the case of the Feynman rules in the particle physics context, 
but, unfortunately, is not very well known by the cosmology community. In the present work, we extend such 
' ' an approach in order to include not only scalar field perturbations as the generators of £, but also vector 

field perturbations. The purpose is twofold: first, we would like the diagrammatic approach (which we would 
call the Feynman-like rules) to become widespread among the cosmology community; second, we intend to 
give an easy tool to formulate any correlator of £ for those cases that involve vector field perturbations and 
that, therefore, may generate prolonged stages of anisotropic expansion and/or important levels of statistical 
anisotropy. Indeed, the usual way of formulating such correlators, using the Wick's theorem, may become 
very clutter and time-consuming. 



1 Introduction 

Our Universe exhibits departures from the exact isotropy and homogeneity; we may see it from the distribution of 
temperatures in the cosmic microwave background radiation (CMB) [1] and from the distribution of matter that 
gives rise to the large-scale structure [5] • The origin of such departures is a marvelous phenomenon that involves 
two outstanding mechanisms in particle cosmology [3J: the production of virtual particles via the Heisenberg 
uncertainty principle, and the primordial accelerated expansion of the Universe that defines a dynamical particle 
horizon. The joint work of these two mechanisms leads to the production of real particles (the Hawking radiation 
phenomenon) and the classicalisation of the field perturbations that live in the spacetime which describes our 
Universe. The nature of such field perturbations may be scalar, vectorial, spinorial, or even, it may correspond 



*e-mail: cesar . valenzuelaScorreounivalle . edu . co 
T e-mail: yeinzon. rodriguezOuan. edu. co 
'-e-mail: juanpbeltran@uan.edu.co 



to p-form a^ 1 ! . In addition, their contribution to the departure from the exact isotropy and homogeneity (pa- 
rameterized by the primordial curvature perturbation C El IS]) m ay be operative during inflation (the standard 
inflationary mechanism |T01 [TTJ [T3J H31 [El E3 ) > at the end of inflation (the inhomogeneous reheating mechanism 
j!6j). or after inflation (the curvaton mechanism [T7J HH [H] ) • Depending on the nature of the field perturbations, 
on the number of fields involved, and on the action describing the theory of gravity and the matter content of 
the Universe, the statistical properties of the distribution of temperatures in the CMB and of the distribution of 
matter that form the large-scale structure may exhibit different features: once the variance has been fixed, the 
power spectrum of £ may become scale-dependent, either blue-tilted or red-tilted, or even may become dependent 
on the direction of the wavevector (signaling violation of the rotational invariance in the two-point correlator of 
£, a property which is called statistical anisotropy [4J [20]). The higher-order correlators of £ are also affected: the 
odd-point correlators may become different to zero (signaling departures from a gaussian distribution [§]) and the 
even-point correlators may become different to products of two-point correlators (signaling also non-gaussianity 
[3]). Moreover, the three-point correlator may depend on the directions of the three wavevectors, something 
which has been called anisotropic non-gaussianity |2I1 I22| . Talking a bit more about the statistical anisotropy, it 
is generated by either a primordial anisotropic expansion, or a non-scalar nature of the fields involved, or both. 

Calcula ting the n-point correlators of £ from a well-defined action is crucial in order to compare theory and 
observation^ 2 !. The right framework to propagate the statistical properties of the field perturbations to the 
statistical properties of £ is the cosmological perturbation theory (CPT) [H [23j [24l [25] ; however, this normally 
involves lengthy calculations, even more if the nature of the fields is not scalar [3S] [371 (Ml US]- In addition, 
interesting phenomena such as non-gaussianity may only be accessible if the CPT is taken to second or higher 
orders [9] [30] [31]. Despite of this, the CPT is valid throughout all scales, leaving no room for discrepancies 
attributed to not considered subhorizon phenomena. A different approach is via the SN formalism [32 , 33 , 34 , 35 , 
where £ is identified with the perturbation in the number of e-folds of expansion from an initial time in a flat 
slicing to a final time in a uniform energy density slicing (the threading must be comoving). The SN formalism 
gives an expression for £ which is valid to all orders in CPT; however, it is only valid for superhorizon scales (in 
absolut contrast with CPT) [9] . Of course, extracting the statistical properties of the distribution of £ in the SN 
formalism requires also to do some "perturbation theory" : to expand £ = SN in a Taylor series and to cut it out 
at the desired order [31] [36] . The classicalisation is given at superhorizon scales [37] , as well as the conservation 
of £ is if the adiabatic pressure condition is hold [9] [35] ; that is why we do not worry much about the subhorizon 
scales and still we obtain very precise results at a very cheap computational cost. 

The way of calculating n-point correlators of £ in the SN formalism is very similar to the way of calculating 
scattering amplitudes in the cannonical formulation of Quantum Field Theory [38]: the Wick's theorem is 
essential [39], the calculation is very direct but not so intuitive, and anyway it may become very clutter and 
time-consuming. Feynman was very clever, even for minor aspects such as a void ing not so intuitive procedures 
when calculating things; he developed a series of diagrammatic rules [40] I41)# 3 I that can reproduce every piece 
of the calculation of scattering amplitudes, very easy to apply, very easy to remember, and very illuminating at 
trying to visualize what is actually happening in a quantum process. Feynman rules became widespread among 
the particle physicists community and it is nowadays the standard language to interpret quantum processes and 
calculate scattering amplitudes [43] . In view of the above, a few years ago a diagrammatic approach to calculate 
n-point correlators of £ was developed including only scalar fields as the generators of £ }44)# 4 I . Such "Feynman- 
like" rules share the same advantages of its counterpart in Quantum Field Theory, making very easy and intuitive 
to start any n-point correlator calculation. External lines in this case are identified as Fourier modes of £ and 
internal lines, the propagators, are identified as the power spectra of the scalar field perturbations. If the field 
perturbations are non-gaussian, new diagrams appear which do not have a correspondence in the particle physics 
Feynman rules: internal lines that split like branches of a tree; they correspond to connected n-point correlators 



* 1 The particle production process has been well studied in the scalar field case [3]. Some nice studies in the vector field case may 
be found in Refs. dEHIDE]. 

# 2 At the end, the theory cannot predict the exact outcome of an experiment in a member of the ensemble since the quantum 
nature only allows us to predict probabilities in the ensemble. The ergodic theorem is crucial in this respect [8] and, therefore, the 
assumption of statistical homogeneity (translational invariance of the n-point correlators) must be maintained. 

# 3 Actually, this technique was first described by Feynman at the Poconos Conference in 1948, and published in 1951 1421 . 

^ 4 A closely related diagrammatic formalism was developed some time ago in order to study large-scale structure formation via 
gravitational instability; see for instance Ref. 1451 and references therein. See also Ref. 1461 for a diagrammatic treatment of 
perturbative calculations in the "in- in" formalism and Refs. 1471 1481 for a diagrammatic approach to information field theory applied 
to the reconstructions of non-linear signals. 



of the field perturbations. The issue of the vertex renormalization was also treated, showing in an elegant way 
how to absorb the diagrams with dressed vertices int o ju st one diagram with undressed, but redefined, vertices, 
in complete analogy with the Feynman rules case [43J# 5 I. 

In this paper, we want to draw the attention of the reader into two significant aspects of the history of the 
Feynman-like diagrammatic approach in cosmology. First, although very powerful as it is in the particle physics 
context, not many cosmologists are aware of the existence or the power of the Feynman-like rules, and still 
most of the calculations based in the 5N formalism are performed directly, invoking the Wick's theorem. Of 
course, the arrival to the relevant integrals is lengthy, the intuitive connection between an important feature in 
the calculation and the physics behind it is almost lost, and the scene quickly becomes very clutter. Second, 
vector field perturbations were not considered in the original formulation of the Feynman-like rules. Vector fields 
as possible generators of £ have attracted the attention of several scientists as they may generate significant 
statistical anisotropy, still in agreement with present observations [HJ El Ell E2] i via the anisotropic expansion 
they may generate gS W\ 121 US E3] and also via its own vector nature HEJEHHEIEIEIEIIEE]- Indeed > 
statistical anisotropy surely will become discriminator of models for the generation of £ taking into account the 
forthcoming increment in the precision of observations starting with the PLANCK satellite data [H [59] . The 
first purpose of this paper is to try to disseminate the diagrammatic approach to cosmologists as much as possible; 
time will tell if this objective was reached. The second purpose of this paper is to extend the Feynman-like rules 
to the case where the contributions to £ are due to several scalar and vector fields, producing the latter, in 
general, anisotropic expansion, statistical anisotropy, and anisotropic non-gaussianity [2U EOl Ell E21 (Ml (Ml E3 . 
We try to be as general as possible, allowing for an arbitrary number of fields, taking into account the transverse 
as well as the longitudinal polarizations of the vector fields (i.e. allowing for mass terms), keeping in mind that 
the action may not be parity-invariant, considering, in general, anisotropic expansion, and assuming that the 



probability distribution functions for the scalar field perturbations, as well as for the vector field perturbations!^ 
may in general be non-gaussian. It is worth emphasizing that, in the context of nowadays cosmology, it is very 
important to have an efficient method for the calculation of loop corrections; this is because high precision 
cosmological measurements could be able to access observational signatures encoded in the loop diagrams if they 
turn to be non- negligible [66], so it is crucial to go beyond tree level terms. Moreover, from the theoretical point 
of view, there exist scenarios in which loop contributions are larger than the tree-level terms and, therefore, 
are able to produce large non-gaussianities [60l EU [67l El ESI EH); in those scenarios, it is essential to have an 
efficient method to evaluate such contributions. 

The layout of the paper is as follows. In Section 2 the n-point correlators are introduced; the meaning of 
statistical homogeneity, statistical isotropy, and gaussianity is clearly established. In Section 3 the 8N formalism 
is introduced, allowing for the possibility of anisotropic expansion and, therefore, for the existence of vector 
fields as generators of £; the two-point correlator of C, is calculated using the non-diagrammatic technique (i.e. 
by employing the Wick's theorem) up to the one-loop level. Section 4 is devoted to the presentation of the 
Feynman-like rules in the limit where the field perturbations are gaussian; the calculation of the two- and three- 
point correlator of C, up to the one-loop level, and four-point correlator up to the tree level, are performed 
employing the diagrammatic approach and compared with the findings in Section 3. Section 5 introduces the 
rules in the non-gaussian case; again, some examples are considered and compared with the non-diagrammatic 
approach. In Section 6 the vertex renormalisation is studied for the gaussian and non-gaussian cases. The 
conclusions of this work are presented in Section 7. 

2 Statistical homogeneity, statistical isotropy, and gaussianity 

Quantum mechanics only allows us to predict probabilities of different outcomes after an experiment in an 
ensemble of systems, in contrast to classical mechanics which does allow us to predict the exact outcome after an 
experiment in just one element of the ensemble. Since the underlying physical mechanism in the generation of 
cosmological perturbations is of quantum nature, the cosmologists are more interested in studying the statistical 
properties of a perturbation map, say the CMB map or the galaxy distribution map. The way of doing this is 
via the n-point correlators of the perturbations in real space. Let's define a scalar cosmological perturbation 



# 5 The standard vertex corresponds to a m-order derivative (being m the number of internal lines ending in that vertex) of the 
unperturbed number of e- folds N with respect to the fields involved 1441 . 
# 6 Actually, the probability distribution functions for the scalar perturbations that multiply the respective polarization vectors. 



/3(x) in real space and make a Fourier integral expansion 

/rl 3 k 
^3 e^/3(k) , (1) 

where /3(k) is the Fourier mode function of /3(x). The n-point correlators of /?(x) are aver ages over the ensemble 
of the products /3(xi)/3(x 2 ).../3(x„) where Xi, X2, ..., x n represent different points in space 



#T 



(/3(x 1 )/3(x 2 ).../3(x„)) ee /^L^%...|!%e^ kl - +k — +- +k "- x ")(/3(k 1 )/?(k 2 ).../3(k„)) . (2) 



(27r) 3 (27r) 3 '"(27r) 3 

Thus, the correlation functions in real space may be studied via the correlation functions in momentum space. 
Let's see now the meaning of statistical homogeneity, statistical isotropy, and gaussianity. 

2.1 Statistical homogeneity 

Of course the perturbation map is not homogeneous (i.e., it is not invariant under spatial translations), but 
it may be that the probability distribution function governing /3(x) is, which is called statistical homogeneity 
[51 [i?l[2"D]. This means that the n-point correlators in real space arc invariant under translations in space, i.e. 

G0(xi + d)/?(x 2 + d).../3(x„ + d)) = (/3(x 1 )/3(x 2 ).../3(x„)) , (3) 

where d is some vector in real space establishing the amount of spatial translation. The only way of achieving 
this, in view of Eq. ([5]), is expressing the argument in the exponential function inside the integral as the addition 
of several terms of the form /(xj — Xj), which in turn is possible (but it is not the only possibility) if the n-point 
correlators in momentum space are proportional to a Dirac delta function: 

(0(ki)0(k 2 )...0(k„)) = (2 7 r) 3 ( 5 3 (k 12 ... n )A/ /3 (k 1 ,k 2 ,...,k„). (4) 

In the previous expression, ki 2 ...„ means ki + k 2 + ... + k„, and the function Af / g(ki, k 2 , ...,k n ) is called the 
(?i — l)-spectrum. Statistical homogeneity is absolutely necessary as an hypothesis of the ergodic theorem [8j; 
otherwise, although we do have a theoretical framework to do calculations and although we do have a significant 
amount of observational data with unprecedented precision, we could not compare one with the other. 

2.2 Statistical isotropy 

Once statistical homogeneity has been secured, in the form of Eq. Q, we ask about the invariance under 
spatial rotations (i.e. isotropy). Of course again, the perturbation map is not isotropic, but it may be that the 
probability distribution function governing /3(x) is, which is called statistical isotropy [8j[9][20]. This means that 
the n-point correlators in real space are invariant under rotations in space, i.e. 

( i 9(xi) j 0(x a )... i 9(x n )) = (/3(x 1 )/3(x 2 ).../3(x„)) , (5) 

where x, = 1Z x^, 1Z being a rotation operator. To satisfy the above requirement, the (n — l)-spectrum must 
satisfy the condition 

M^ki, k 2 , ..., k„) = M^(k!,k 2 , ..., k„) , (6) 

where the tildes over the momenta represent as well a spatial rotation, parameterized by TZ, in momentum space. 
This condition has more explicit consequences in the spectrum (1-spectrum) and the bispectrum (2-spectrum): 

M /3 (k 1 ,k 2 ) = P / ,(ki > k 2 )=P / j(fc) > (7) 

M^(ki,k 2l k 3 ) = ^(ki > k2,k 3 )=B^(fci,fc2 J A3), (8) 

where in the first line k = |ki| = |k 2 |, and in the second line ki = |kj|. Starting from the trispectrum (3- 
spectrum), the condition in Eq. (j6]) about statistical isotropy in all the higher-order (n — l)-spectra cannot be 
reduced to similar conditions to the ones in Eqs. and (|8]l. so that the minimal way of parameterizing the 



# 7 The ensemble average inside the integral is over the Fourier mode functions only since they are the stochastic variables. 



(n — l)-spectra (with n > 4) will always be in terms of all the n wavevectors. The scalar nature of /3(x) is very 
important since, if it were a vector or a tensor, there would not be a way to make the n-point correlators in real 
space invariant under spatial rotations. In those cases, we relax the meaning of statistical isotropy and establish 
that it is present if the (n— l)-spectra of the scalar pertubations that multiply the respective polarization vectors 
or tensors satisfy Eq. (|6|). 

2.3 Gaussianity 

Gaussianity may be defined either via the perturbations /3(x) in real space or via the perturbations /?(k) in 
momentum space. Let's talk first about the latter and we will come back later on the former. We say that the 
probability distr ibution function governing /3(k) is gaussian if for different wavevectors the perturbations are 
uncorrelatedl #8 l 

(PQcJpfa)) = (2^) 3 <5 3 (k 12 )P^(k 1 ) , (9) 

and if the n-point correlators, with n being odd, are zero, while those with n being even are e qual to the sum 



over all ways of pairing /?(k)s with each other of a product of the 2-point correlators of the paira^ 9 ! [El [H] : 



(/3(k 1 )/3(k 2 )/3(k 3 )) = 0, (10) 
(/3(k 1 )/3(k 2 )/3(k 3 )/3(k 4 )) = (/3(k 1 )/3(k 2 ))(/3(k 3 )/3(k 4 )) + (/3(k 1 )/3(k 3 ))(/3(k 2 )/3(k 4 )) + 

(/3(k 1 )/3(k 4 ))(/3(k 2 )/3(k 3 )) (11) 

= (2 7 r 6 )(5 3 (k 12 )c$ 3 (k3 4 )P /3 (k 1 )P / g(k 3 )+ two permutations, (12) 

and so on. Eqs. (|9|) and (fT2"j) clearly show that statistical homogeneity is a necessary but not sufficient condition 
if the probability distribution function governing /3(k) is to be gaussian [HE]. Now, coming back to the pertur- 
bations /3(x) in real space, it is possible to show that the gaussianity condition expressed above implies that the 
probability distribution function P(/3(x)) is given by 

P(/3(x)) = ] .. e-^WW> , (13) 

VM^ ( x )) 

which is the usual definition of a gaussian probability distribution function. However, requiring Eqs. (jlpp . (jlll) . 
and so on for /3(x), instead of /3(k), also results in Eq. (|13|) without requiring gaussianity in the mode functions 
/?(k). In other words, gaussianity in /3(k) is a sufficient but not necessary condition for gaussianity in /3(x); 
indeed, via the central limit theorem [71] [72], /3(x) may become gaussian just by being expressed as a sum of 
uncorrclated quantities (/3(k)) even if they are not gaussian. 

When the perturbations /3(k) arc non-gaussian, but still assuming statistical homogeneity, the n-point cor- 
relators in momentum space, with n > 3, are expressed in terms of "connected" n-point correlators (identified 
with a subscript c) which establish the departure from the gaussianity condition: 

(/3(ki)/3(k 2 )) = (2ir) 3 S 3 (k 12 )P f! (k 1 ), (14) 

(/3(k 1 )/3(k 2 )/3(k 3 )) = (/3(k 1 )/3(k 2 )/3(k 3 )) c (15) 

= (2^) 3 <5 3 (k 123 )P^(k 1 ,k 2 ,k 3 ), (16) 
(/3(k 1 )/3(k 2 )/3(k 3 )/3(k 4 )) = (/3(k 1 )/3(k 2 )/3(k 3 )/3(k 4 )) c + (/3(k 1 )/3(k 2 ))(/3(k 3 )/3(k 4 )) + 

(/3(k 1 )/3(k 3 ))(/3(k 2 )/3(k 4 )) + (/3(k 1 )/3(k 4 ))(/3(k 2 )/3(k 3 )) (17) 
= (2^) 3 ( 5 3 (k 1234 )T /3 (k 1 ,k 2 ,k 3 ,k 4 ) + 

(2 7 r 6 )(5 3 (ki 2 ) ( 5 3 (k 34 )P / 3(ki)P ; 3(k 3 ) +two permutations, (18) 

and so on. In the above, P^(ki,k 2 ,k 3 ) and T^ki, k 2 , k 3 ,k 4 ) are called the connected bispectrum and trispec- 
trum of j3. 



# 8 For equal wavevectors, the reality condition on /3(x) applies and, therefore, there is self-correlation. 

* 9 The sum over pairings does not distinguish those which interchange wavectors in a pair, or which merely interchange pairs. 



3 The primordial curvature perturbation and the 5N formalism 

The SN formalism 32, 33, 34, 35 provides a powerful method to evaluate the primordial curvature perturbation 
C(x, t) in terms of the perturbations of the fields at the the time of horizon crossing t* (corresponding to a flat 
slicing), and the derivatives of the unperturbed number of e-foldings 7V(t,t*) = J, H(t')dt' with respect to the 
unperturbed fields evaluated at £*. 

According to this formalism, once the separate universe approach has been invoked (see also Ref. |73j). and 
a comoving threading has been established, the value of (, in a uniform energy density hypersurface at the final 
time t is given by the perturbation in the time integral of the local volume expansion rate 9 along a curve starting 
at an initial flat hypersurface at the time U: 

C(x, t) = 5N(x, t, U) - (5N(x, t, U)) . (19) 

Here, the bracket notation means a ensemble average (which corresponds to a spatial average is there is statistical 
homogeneity). In many inflationary scenarios, the number TV of e-foldings depends only on the values of the 
fields at horizon crossing so we can write the curvature perturbation as an expansion in the perturbations of the 
fields at this time [36]. Supposing that inflation is driven by a single scalar field <fi(x,t) which can be decomposed 
as <f> = 4>o + 5(f> where S(f> is the field perturbation at the initial time £, = t* , the <5TV formula acquires the form: 

,. . dN(t,U). ±l , 1 <9 2 TV(i,i*) rj2 , , 1 d 3 N(t,U). l3 , , . . , nn , 

cm = -iJL-Lsfat.) + ^^r^ 2 ^^) + 3? V(tr^ (x '^ } + (spatial average) ■ (20) 

By doing this expansion, we are neglecting any additional dependence on <fi after horizon crossing which implies 
that we assume that the field perturbations are strongly suppressed on large scales; this assumption is valid 
during a slow-roll regime but, in general, it is valid for any light fields present during inflation. To simplify the 
writing of the expressions, in the following we omit the dependence of the field perturbations on the initial time 
t* ; we also omit the spatial average in the 5 TV formula. 

Due to the versatility of the method, it has been implemented in many different scenarios; for instance, if we 
assume that the primordial curvature perturbation is generated by n scalar fields, the <5TV formalism gives us the 
following expansion for (: 

C(x, t) = SNfa (x), . . . , 0„(x), t) = NrfMx) + ^NuSMxWA*) + ^ JJK ^(x)^ J (x)^ if (x) + . . . , (21) 

where Nj = dN/d(f>i, etc. / = 1, 2, • • • , n and we understand that repeated indices are summed over. In a recent 
work this formalism was extended to include also vector field perturbations [I] (see also Ref. [12]), and it was 
shown that in the simplest case where £ is generated by a single vector field and a single scalar, the curvature 
perturbation can be calculated by means of the following expression: 

C(x, t) = 6N(cf>(x), Ai(x), t) = N+dt + NidAi + ^N^(Scj)) 2 + N^S^bSA, + ^NySAiSAj + ..., (22) 

where 

<9TV _ <9 2 TV d 2 TV , . 

TV,= — , N„ = w , TV, J = ^^, etc., (23) 

are the derivatives with respect to the scalar <j> and the spatial components of the vector field A. It is very 
important to warning the reader that the extension of the 5N formalism to the vector field perturbations case 
requires relaxing one assumption in the separate universe approach: the isotropic expansion. The consequence 
of this is to have a tensor perturbation hij which is time-dependent [4]; however, the identification of £ with SN 
remains the same no matter that hij is time-dependent or not. If for some reason the expansion is isotropic, 
the perturbations in the scalar and vector fields may be statistically isotropic; however, £ would not be in 
general statistically isotropic due to the vector nature of the fields involved (except, of course, for <T>). On the 
contrary, if the expansion is anisotropic, all the perturbations would be automatically statistically anisotropic 
and, therefore, ( would be too; this would affe ct al so the TV derivatives since the background would not be of 
the Friedmann-Robcrtson- Walker type anymore ^ 10 ! . 

# 10 N depends on the values of the relevant fields at the initial time t; and on the total (and uniform) energy density p at the 
final time t; there is no explicit dependence on the position. The TV derivatives are evaluated in the background so there is neither 
explicit nor implicit dependence on the position unless the background metric is inhomogeneous or some very specific configuration 
of inhomogeneous background fields renders the background metric homogeneous; however, in the former case, there would be 
statistical inhomogeneity, making impossible to compare theory and observations. 



Now, we move to the more general case in which the inflationary dynamics is driven by multiple scalar and 
multiple vector fields. We will consider n scalar fields and m vector fields. To deal with the contributions coming 
from the different fields involved, we introduce the notation: 

5$ A = {60j , 5A1) . (24) 

The index A is separated in two sets, a set of indices I labelling the scalar fields which runs from 1 to n and 
another set of indices a labelling vecto r fie lds which runs from 1 to m. The index i specifies the component of 



any vector field and it runs from 1 to 3 W I . We use the index a as a supra- index and the index i as a sub- index. 



Accordingly, the derivatives of N with respect to the fields are separated as follows: 

N A = {N!,N?}, (25) 

Nab = {Nu , Nj , N$} t (26) 

and so on. For instance, in the notation above, we represent the mixed second derivative with respect to 4>i and 
A b , as Njj = d 2 N/d4>idA t '. In terms of the notation in Eq. (|24[) . the curvature perturbation for the multi-scalar 
and multi- vector field case is written as 

C(x,t) = <W($5(x),i) 

= N A 8$ A + -N AB S^ A S^> B + j l N ABC 6$ A 5$ B 5$ + -N abcd 8<5> a 5<5>b^c^d + ■ ■ ■ ■ (27) 

As we can see from the equations above, since £ is a scries in the field perturbations 5& A , then its correlation 
functions will be expressed as perturbative series in the correlation functions of 5& A . To this end, we first write 



the mode function associated to ( in terms of the mode functions associated to the field perturbation a^ 12 



1 f r/ 3 A-i 

C(k,t) = N A 5$ A Qt) + -N AS j — ^<5$ 1 (k-k 1 )^ s (k 1 ) + 

T f s$ A (k - ki - k2)<y*s(kx)^ e (k 2 ) 



i AT f d 3 k lC i 3 k 2 

W^ ASC J (2tt)6 ' 



^N AS (2Tr) 3 6 3 (k)(5$ A (x)S<S> s (x)) - -N ASC (2ir) ;i 5 3 (k)(5<S> A (x)6<S> s (x)6<S>c(x)) - ... , (28) 

where in the last line the mode function of the spatial average in the 5N formula is shown. Now, we may multiply 
n mode functions of £ and perform the average over the ensemble in order to obtain the n-point correlator. As 
an example, we will obtain the 2-point correlator of £ by truncating the series up to third order and taking into 



* 11 This is because, as inflation homogenizes the vector fields (diA^ = 0), we expect the temporal components of them to vanish 
[5]. Besides, if the vector fields are massless, we can set the temporal components to zero by a gauge choice. 

* 12 Here we take into account Footnote |#10| even for those cases where the background metric is anisotropic, but still homogeneous, 
the N derivatives evaluated in the background do not depend either explicitly nor implicitly on the position. 



account terms that will only be represented by tree level or one-loop diagram s^ 13 
(C(ki)C(k 2 )) = N A N B (8$ A (k 1 )8$ B (k 2 )) + 



\n a n bc 



^NaNscd 



d 3 k 3 



^(<5$ A (k 1 )5$ i5 (k 2 - k 3 )<5$ c (k 3 )) + 

(<J* yi (ki)<y$ fl (k 2 - kg - k 4 )5$ c (k3)^ s (k4)) + 



(27T) 

d 3 k 3 d 3 k 4 



(27T) 



1 



■7 N abN S d 



d 3 k 3 
(2tt) 3 
d 3 k 3 d 3 k 4 



(27T) 



(<5$ A (k 2 )<5$ i j(k 1 - k 3 )5$ c? (k3)> + 

(<J$x(ki - k 3 )<5$ i 3(k3)5$ ( 5(k2 - k 4 )<5$^(k 4 )) 



I /" (V 3 A-Q 

- A N AS N CD {2n?5 3 {\^) J -^(8§ A {k x -k3)^g(k3))<^e(x)^s(x)> + 

j3r, j3i. 

-5- r i( ( j$ y i(k 2 )<y$ fi (ki - k 3 - k4)<J$c(k3)^s(k4)> - 



(27r)e 



qN a N B0S 

^AB^cB(2^) 3 ^ 3 (ki) 
1 



rf 3 fc 3 
(2tt)* 



(<5$ A (k 2 - k 3 )5$ g (k 3 ))(<5$ c ( X )5$ s (x)) 



iV^7V Cfl (2^) b ^(k 1 ) ( 5 3 (k 2 )( ( 5$^(x)^ s (x))(^ c? (x)5$ z5 (x)) 



(29) 



The following step is to apply the Wick's theorem |39j , assuming that the held perturbations may not be gaussian 
and, again, avoiding terms that would be represented by two- or higher-order loop diagrams: 



(C(kx)C(k 2 )) 



N A N B {8$ A {k x )8$ B {k 2 )) 



\n A Nbc 



yN A N S0 D 



d 3 k 3 



(5^(kx)5$ s (k 2 -k 3 )5$c(k 3 ))c 



(2tt) 3 

d 3 k 3 d 3 k 4 ' 



{2ttY 



(<5^(ki)<5$ s (k 2 - k 3 - k A ))(S$>c(k 3 )S<S> B {k 4 )) 



(<5^(k 1 )^ e (k3))(^s(k 2 - k 3 - k 4 )^(k 4 )) 
< < J$ ji (k 1 )<J$ fl (k4))(5$ s (k 2 - ka - k 4 )^ (5 (k3)) 



^N A N BG 



d 3 k 3 



1 



-NabNcd 



(2tt) 3 
d 3 k 3 d 3 k 4 



(5$^(k 2 )5$ S (ki-k 3 )5$c(k3))c 



4"'""'°"./ (2tt) 6 

(<J^(k! - k 3 )<5$c(k 2 - k 4 ))(<5$B(k 3 )^ s (k 4 )) 

(<J#x(ki ~ k 3 )8$ B (k 4 ))(8<S> B (k 3 )8$ c (k 2 - k 4 )) 



(5^(ki - k 3 ) < 5$B(k 3 ))(J$ c (k 2 - k 4 )<S$ B (k 4 )) 



\N AS N CB (27r) 3 8 3 (k 2 )(8<S> c ,(x)S<S> B ( X )) 
1 f d 3 k 3 d 3 k A f 

yN A N B0B J -j^y e 



d 3 k 3 

(2ny 



(^ 1 (kx-k 3 )^ 5 (k 3 )) 



( ( 5$^(k 2 )J$ i5 (k 1 - k 3 - k 4 ))<5$ c (k 3 )<5$ fl (k 4 )) 



(<5$^(k 2 )^<5(k3))(^B(ki - ka - k 4 )<r$ s (k 4 )) 

(,5$ j 4(k 2 )5$ B (k 4 ))(5$ s (k 1 - ka - k 4 )<5$ c (k 3 )) 

1 ( d 3 ki 

-^ AB N 5B (2^ 3 8 3 {k 1 )^ 5 ^)8^ B ^)) J j^± 

7^s^CD(2^) 6 ^(ki)<5 3 (k 2 )(5$^(x)5$ s (x))(5$ ( 5(x)5$ 5 (x)). 



(8$ A (k 2 -k 3 )8<!> B (k 3 )) 



(30) 



# 13 Except for four terms which are formally represented by two-loop disconnected diagrams and that, as we will see, will cancel 
each other out. 



We now employ the definition in Eq. (|14[) for the power spectrum of £, and the expressions 

(,5<Mki)M> s (k 2 )) ee (2^) 3 <5(k 12 )n Ai5 (k 1 ), (31) 

(^^(k 1 ) ( 5$ s (k 2 )^ C7 (k3)) c = (27r) 3 6(k 123 )B ASC (k 1 ,k 2 M), (32) 

(5$ >i (k 1 )^fl(k 2 )«$c(k3)**s(k4)>c = {2^) 3 5{k l23i )T ASCD {k 1 ,k 2 MM), (33) 

and so on; these expressions define the power spectra H^B °f the field perturbations and the higher order 

functions BabCi ^Abcd^ an d so on, that parameterize the intrinsic non-gaussianity in the field perturbations. 

Thus, calculating some trivial integrals that involve Dirac delta functions we arrive to 



(2 7 r) d ^(k 1 +k 2 )P c (k 1 ) =(27r) d a*(ki 

1 



k 2 

^ a N BC d^ab(^i 
(2n) 6 S 3 (k 1 )5 3 (k 2 ) 
(2 7 r) 3 (5 3 (k 1 +k 



NaNb^ab(^i) + \n a Nsc 



d 3 k 3 
(2nf 



■ B ABc(ki,k2 -k 3 ,k 3 ) + 



d 3 k 3 



n cs (k 3 ) 

d 3 k 3 d 3 k 4 



(2tt) 3 " 
NmN, 



N A N 3 c 



d 3 k 3 

(2tt) 3 



Babc^^i ~k 3 ,k 3 



"AB^CD 



\nabNcd 

1 



(2tt)« 

d 3 k 3 



n^(k 1 -k 3 )n c?5 (k 2 -k 4 ) 



+ 



(2tt)* 



n^(ki-k 3 )n^(k 3 ) 



(27r) b ,5 3 (k 1 ) ( 5 3 (k 2 ) ^N AB N OD (8$c(x)5$ D (x)) 

d 3 k a 



(27r) 3 ^(k 1+ k 2 ; 



-N A N S c B U AS (k 2 ) 



(2nf 



n 



(2 7 r) & ,5 3 (k 1 ) ( 5 3 (k 2 ) -N Ab Nc5(S$c(x-)6*d(x-)) 



-N AS N CD (S<S> c (x)5<S> D (x)) (6$ A (x)6$ s (x)) 



d 3 k 3 
{2-kY 

CD( k 3) 

d 3 k 3 

(2tt) 3 



n As(ki -k 3 



n^(k 2 -k 3 )- 



(34) 

The expression in the third line cancels out exactly the expression in the fifth line; this is because, according 
to Eq. ©, 

/d 3 k 
(^J3 n AB(k). (35) 

The same argument is used to demonstrate that the terms in the last two lines cancel each other out. These 
four terms are the ones that formally are represented by two-loop disconnected diagrams. As we see, having 
defined £ as SN — (SN) was key to make this kind of diagrams disappear. Concluding, we may say that the 
power spectrum of £, up to the one-loop level, is given by 

d 3 k 3 



^(ki 



NaNbUabQh) + \n a n B c 



(2tt) : 



f- B Asc(ki,-ki -k 3 ,k 3 ) 



\ Ni N 
1 



BCD^AS^l) I T^3 U CD( k 3) + ^ N A N BC J 7^3 B ABC (-*■!, k i - k 3 , k 3 ) 



-N Ab Ncd 



^U AC (k 1 - k 3 )U SB (k 3 ) + In a Nscd^ab(-^i) 



d 3 k 3 



The reader may notice that the obtaining of the previous expression for the power spectrum of £ is quite 
lenghty and clutter, as well as direct but not so intuitive. As it may be expected, the higher order correlators 
of £ demand even more time and are definitely very clutter leading to undesirable mistakes if extreme care in 
the manipulation of formulas is not taken; all of this of course means that this methodology of calculation is not 
in any respect efficient. The same situation happened in the cannonical formulation of Quantum Field Theory 
until Feynman presented his diagrammatic rules [40) [41] [42l [43] . Following the same idea, the aim of this paper 
is to present a simple diagrammatic approach to calculate all tree level and loop contributions to the correlation 
functions of £ when vector and scalar fields arc present. We start in Section @] with the gaussian case and we 
treat the non-gaussian case in Section [5] We use the tools and the notations that we have introduced in this 
section throughout the rest of the paper. 



4 Fourier space diagrams for gaussian scalar and vector field pertur- 
bations 

As we discussed in the previous section, the disconnected part of the n-point correlators of C exactly vanishes 
due to the substraction of the ensemble average of SN. We keep this in mind while writing down the rules for 
drawing the diagrams related to the n-point correlator functions of C to the desired order in perturbation theory. 
In this section we present the rules to calculate the n-point correlation functions to p-order (0(P P )) when £ has 
contributions from scalar and vector field perturbations. Following this terminology, the tree level approximation 
of the n-point correlation function has p = n — 1 propagators, the 1-loop approximation has p = n, and so on. 

1. Draw all distinct diagrams with n-external lines and p propagators. External lines are represented by solid 
lines while internal lines are represented by dashed lines. Every vertex connect an external line to at least 
1 propagator. 

2. Label the external lines with incoming momenta k^ for i = 1, ...,n and label the propagator with internal 
momenta pk for k = 1, ...,p. Label each end of each propagator with a field index: A, B, ■ ■ ■ ,C. See Fig. 
ffl 

3. Assign the factor Il^^k) to each propagator in the diagram where k is the momentum associated to the 
propagator which must flow from the end with the label A to the end with the label B (see Fig. [TJ. The 
direction of the momentum flow in the diagram is relevant when the action is not invariant under the parity 
transformation. If the momentum flows from the end with the label B to the end with the label A, the 
assigned propagator must be Ilg^k) = Hjl§(— k). 

4. Assign a factor N^ £ .£ (27r) 3 <5(ki — pi... t ) to each vertex, a i-vertex, where the number t of derivatives of 
N is the number of propagators attached to this vertex. We use the convention that incoming momentum 
is positive. The Dirac delta function ensures that the momentum is conserved. See Fig. [5J 

5. Integrate over the propagator momenta ,A 3 J d 3 pi. The first n — 1 integrals can be done immediately 
using the Dirac delta functions but any further integral in general cannot be performed analytically. This 
is the case when there are integrals corresponding to loop corrections. 

6. Divide by the appropriate numerical factor: 

• l\ if there are I propagators attached to the same vertices at both ends. 

• 2 Z Z! if there are I propagators with both ends attached to the same vertex. When the two ends of a 
propagator attach the same vertex, we say that the propagator is dressing the vertex. See Fig. [3] 

7. Add all permutations of the diagrams corresponding to all the distinct ways to relabel the k^ attached to 
the external lines. The number of permutations depends on the symmetries of the diagram: a diagram 
which is totally symmetric with respect to external lines has only one term; in contrast, a diagram without 
symmetries with respect to external lines has n! permutations. 






n„(k) J n LI (k) iV<k) V(k) V^ 

->AAA/ > W\^-- > A/\AAA 

a b 



Figure 1: The diagrammatic representation of the propagator: for each pair of indices A, B we draw a thick 
dashed line. In order to correctly write the propagator LI^g(k), the momentum must flow from the end with 
the label A to the end with the label B; this is relevant if the action is not parity-invariant. In this figure, we 
show the decomposition of the general propagator in terms of one scalar-scalar propagator, two scalar-vector 
propagators, and one vector- vector propagator. Thus, for each pair of scalar indices /, J we draw a dashed line, 
for each pair of scalar- vector indices /, b we draw a dashed-curly line, and for each pair of vector indices a, b we 
draw a curly line. For each vector index a we understand that there is implicitly an index i corresponding to the 
components of the vector field. The diagrams in this paper are drawn using JaxoDraw |74[ 175) . 
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(27i) 3 5(k-q 1 ,)N iR (2^) 3 5(k-q 12 )N u (2;r) 3 5(k-q 12 ) N,. b (2ti) 3 5(k-q 12 ) N.. ab 

— » — < > 




J" 



q 2 X 



* t, ' *" 

q\ 




Figure 2: The diagrammatic representation of the 2-vertex N^g- Momentum must be conserved at each vertex: 
the Dirac delta function guarantees it. Analogously to Fig. [TJ the 2-vertex has been decomposed in the scalar- 
scalar, scalar- vector, and vector-vector 2-vertices. 
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%^r* 








y 
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1 




o 






+ - 
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-KfKSKfKr^ 


+ 


— ^C/N/V/V/V; 



Figure 3: A one-loop Feynman-like diagram for power spectrum P^ with gaussian scalar and vector field perturba- 
tions (in the decomposition we have assumed isotropic expansion so that there are no scalar- vector propagators). 
We call it a "dressed vertex" when the two ends of a propagator are attached to the same vertex as the left one 
in this figure. 

In the rules above, we have two basic building blocks. The first one is the Feynman-like propagator n^g(k): 

iWk) = (n 7J (k) , irj(k) , n^(k) , n« 6 (k)} , (37) 

where the II/j(k) are the scalar-scalar field perturbations power spectra defined by 

(ty/(ki)W(k 2 )) = (27r) 3 <5 3 (k 1 - k 2 )n JJ (k 1 ) , (38) 

where we suppose that in general the different scalar fields could be correlated at horizon crossing so that they 
do not necessarily have the same power spectrum Ps<j>', they depend on the wavevector since , in general, there 



will be anisotropic expansion, rendering the field perturbations spectra statistically anisotropic !^ 14 l . On the other 



hand, the n 7 J(k) are the scalar-vector field perturbations power spectra related to the two-point correlators of 



a scalar field <fii and a vector field A b whose origin lies again in the anisotropic expansion 

T b 



(S^ I (k 1 )SA b 1 (k 2 )) = (27rf5 3 (k 1 + k 2 )n '(ki) . (39) 



Here, n 7 ^(k) is given by: 

long 



lC S (k)P long (k), (40) 

where it has been written in terms of the longitudinal component of the power spectra P io b and the parity 



conserving and violating power spectra P + and P_ respectively: 

P±=l(PR±P L b )> (41) 

where Pjj and P L denote the power spectra for the transverse components with right-handed and left-handed 
circular polarizations, and again all of them depend on the wavevector because of the anisotropic expansion. The 



# 14 The particle production mechanism is statistically anisotropic. 

11 



formal definitions of the polarization spectra are given by 

(6^(^)6 A\*(k 2 )) = (2n) 3 5 3 (k 1 - k 2 )P A 6 (k 1 ) , (42) 

where A denotes the different polarizations, L, R, or long, and 

^k)=£e*(k)<L4 A (k), (43) 

A 

e A (k) being the respective polarization vector, and a hat denoting a unit vector. The basis IIj™ n , IIj dd and 

n 7 ° ng is given by 

II!7*(k) = y/2pi , D?f (k) = V2ie jM k k Pi , 1^04= kj, (44) 

where p = IZx and x being the unit vector that goes in the x direction, 1Z being the rotation operator that takes 
z to k, and z being the unit vector that goes in the z direction. Finally, the II".- (k) are the vector-vector field 
perturbations power spectra defined by: 

(M?(ki)M$(k a )) = (2 7 :) 3 S 3 (k 1 + k 2 )n« 6 (k!) , (45) 

where II£ b (k) is given by [4] 

n» b (k) = n?f »(k)pf (k) + in° dd (k)P« 6 (k) + ir^k^jk) , (46) 

^long being the longitudinal component of the power spectra and P+ b and P" b being the parity conserving and 
violating power spectra respectively: 

pab s 1 (pa. ± pa6) > (4?) 

where P^ b and P£ b denote the power spectra for the transverse components with right-handed and left-handed 
circular polarizations; all of them, as in the other cases, depend on the wavevector because of the anisotropic 
expansion. Analogously to the previous case, the formal definitions of the polarization spectra are given by 

(<JAS(ki)^(k3)) = (2 7 r) 3 <5 3 (k 1 - k 2 )Pf (k x ) . (48) 

The basis IIfJ en , n°/ d and n'° ng in this case is given by g] 

n°j cn (k) = 6. tJ - fakj , n° dd (k) = e ^4, n!° ng (k) = hh . (49) 

We represent diagrammatically the full propagator in Fig. Q] 

The second building block of the diagrammatic approach is the t- vertex N% % . The diagrammatic repre- 
sentation of the i-vertex is shown in Fig. [2j There is a factor (27r) 3 i5(ki — P1...4) attached to the vertex. The delta 
function is needed in order to guarantee momentum conservation at the vertex, kj is an external momentum 
(bold solid line), and the pj, are the internal momenta (bold dashed line). 

Examples of the diagrams derived from these rules are discussed in the next section. 

4.1 Examples 

In this section we show how the diagrammatic rules work in a particular case. We will assume that the curvature 
perturbation is generated by a single scalar and a single vector field so that <5<f>^ = (5<j> , 8Aj) and the 5N series 
goes as follows: 

C(x, t) = 8N(</>(x),Ai(x),t) = iV 0( 50 + NiSAi + -N^Sc/)) 2 + N^ i 6(/)SA l + -N.^SA.SAj + (spatial average) . 

(50) 

According to Eq. (|50|) the n-point correlation functions have contributions from the vector field perturbations, 
from the scalar field perturbations and from mixed terms. We will show how the rules work in this case for the 
two-point, the three-point, and the four-point correlators of C- We will work out in detail the case of the 1-loop 
power spectrum and describe schematically the other cases. 
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4.1.1 The tree level and one-loop contributions to the spectrum Pq 

For the power spectrum we have n = 2, hence, the tree level approximation has p = 2 — 1 = 1 propagator. At this 
level we then have only two first order derivatives (1-vertex) and one propagator joining them. Diagrammatically 
the two-point correlator of C is represented by: 



(C(ki)C(k 2 )) t 



+ A/W- 



(51) 



-AAAy + 'WWW- 



Using the rules given at the beginning of this section, and after doing the trivial integration, the tree level term 
to the power spectrum P^ is given by: 

if"(ki) - NjlNbILabQli) = A^(ki) + A^-n^ki) + N^NjIL^-^) + A^H^ki) , (52) 

where P(k) is the scalar-scalar field power spectrum, II^(k) is the scalar- vector field propagator, and ILj(k) is 
the vector-vector field propagator. 

Going now to the one-loop correction, we have n = 2 external lines and p = 2 propagators. We have one 
term with two 2-vertices attached to the two external lines and the two propagators joining the vertices. We 
also have another term with one 1-vertex and one 3- vertex, one propagator joins the vertices and the other one 
dresses the 3-vertex. Assuming isotropic expansion, the one-loop contribution may easily be diagrammatically 
expanded since we avoid scalar-vector propagators; in this case the resulting diagrammatic representation is: 



(C(ki)C(k 2 )> 



1 — loop 




r"v 



■^A/WW^ 



o. 

(53) 



Here we can identify seven different diagrams associated to contributions from scalar, vector and mixed terms. 
Again, using the rules given at the beginning of this section, it is easy to get the one-loop correction to P^ . After 
carrying out one integration (the trivial one) we get: 



P^QH) = ^N AS N GD J ^U AC (k 1+P )U SD (-p) 



d 3 p 

(2tt) 3 

d 3 p 

(2tt) 3 



(27T) 

d 3 p 
^JV|^P(|k! + p\)P(p) + ijVy JVwIIttOd + p)H,-i(-p) + N+iNtjPQk.! + p|)ne(-p) 



\n a Nbcd J ^n^(-ki)n CB ( P ) + \n a n bcd J A^n^(k 1 )n CB ( P ) (54) 



JV^A^P(fci)P(p) + i^^- 00 p( P )n u (-k 1 ) + ijvijv,-^p(p)iiy(ki) 



N^N^P (h^ip) + ^N.N^U^i-k^Uuip) + ^NiN^Uii^Uuip)] ■ (55) 

In the first two lines we have written the total one- loop contribution in terms of the t-vertices N Al ... A and the 
propagators H A § and in the lines below we have expanded them in terms of their scalar-scalar, scalar- vector, 
and vector-vector components. Notice that the first term has a numerical factor 1/2 because, according to rule 

13 



6, there are two propagators attached to the same vertex. The second and third terms also have a factor 1/2 
because there is a propagator dressing a vertex. When we do the expansion of the first term in the scalar and 
vector components the terms with pure-scalar fields and pure- vector fields will retain the numerical factor 1/2, 
however the term associated to mixed contributions (scalar and vector) does not have this numerical factor; this 
is because the propagators are not of the same kind. Similar considerations apply for those terms in the last two 
lines. 

Comparing Eqs. ([36]) . ([52]) and (|54|) we observe that the diagrammatic approach is able to reproduce the 
power spectrum of £ up to the one-loop approximation and assuming that the field perturbations obey a gaussian 
statistics. Nevertheless, the application of the diagrammatic method is very quick, useful, and intuitive, in absolut 
contrast with the traditional method employing the Wick's theorem. Here is where the power of the Feynman-like 
rules relies. 



4.1.2 The tree level and one-loop contributions to the bispectrum B^ 

Assuming £ as in Eq. (|50|) and using the rules given at the beginning of this section to draw the diagrams, the 
tree level term (p = 3 — 1 = 2) to the three-point correlation function of ( is diagrammatically given by: 



(C(k 1 )C(k 2 )C(k 3 )) t 



(56) 



where the diagrammatic expansion has been done assuming isotropic expansion, i.e., avoiding scalar- vector 
propagators. This diagram is the result of joining three external lines with two 1-vertices and one 2- vertex with 
p — 2 propagators. Using the diagrammatic rules we obtain the expression for the tree level bispectrum: 



Bf^kLkij.ka) 



NaNbNcd 



n ic( k i) n 5D( k 2) + 2 perm. 



NlN^[P{k l )P{k 2 ) + 2 perm.] + N t N k N mn n„ n (ki)n fe „(k 2 ) + 2 perm 
N^NiN^j \p(h)U ij (k 2 ) + 5 perm. 



(57) 



The one-loop contribution p — 3 is represented by 



(C(k 1 )c(k 2 )c(k 3 )) 



1 — loop 



• --r 

' / 
' / 

-Vr-' 



/ \ ' 



t \ 



(58) 



which is the result of joining three propagators with three external lines. There are four different ways of joining 
the propagators and the external lines: we have the first term with three 2-vertices, the second term which 
contains one 1-vertex, one 2-vertex and one 3-vertex, the third term with two 1-vertices and one dressed 4- 
vertex, and the fourth term with one 1-vertex, one 2-vertex and one dressed 3-vertex. Again, using the rules and 
performing the trivial integrations (thanks to the Dirac delta function) , the above diagrams lead to the following 
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expression: 



^""(kx.ka.ka) = N AS N 55 N SP J -^ n A5 (-p)n SE (-k 1+ p)n BP (-k 2 -p) 



-N A N BC N BBP 
t:N a N s N CB ef 



\n a N bc N bbp 



d 3 p 



(2tt) 3 
d 3 p 

(2tt 
d 3 p 



(2tt) 3 



n AB( k i) n BB(p) n ci?( k 2 - P) + 5 perm, 
i n5c(ki)n 5 £(k 2 )n^(p) + 2perni. 
n As( k i) n cD( k i2)nB i ?(p) + 5 perm. 



(59) 



We have to expand every term in its scalar-scalar, scalar- vector, and vector-vector components. As an example, 
we expand the first term which has three 2- vertices in the isotropic expansion limit. We will call this term B 222 
because it has three 2-vertices: 



Sf 2 (k l5 k 2 ,k 3 ) 



N aS NcdN BP 



T^p n4 C (-p)iigg(ki + p)n SJ r(k 2 - p) 



/(2^ PWP(|kl+Pl)P(|k2_Pl) 
d 3 p 



N$. 

NijNklNran 



(2tt) 3 
d 3 p 



n ife (-p)rij m (ki + p)n;„(k 2 - p) 
3 {p(p)p(|k 1 + p|)n y -(k 2 -p) + 



P(p)P(|k 2 - p|)IIy (k a + p) + P(|k 2 - p|)P(|k! + p|)n*j(p) 

n^n^n m J -^p |p(p)n ifc (k! + P )n j7 (k 2 - P ) + 

P(|ki + P |)n ifc (p)%(k 2 - P ) + P(|k 2 - pijiiifcCpjnjKki + p; 

The diagrammatic representation of this expansion is shown in Fig. 2J 



(60) 



/ i 
< i 



+ 




+ 



\ 



\ 



Figure 4: The scalar and vector components in the expansion of the B 222 term assuming isotropic expansion. 



4.1.3 The tree level contribution to the trispectrum T^ 

Finally, we draw the diagrams and write down the expression for the tree level term of the trispectrum T^. 
This term has p = 4 — 1 = 3 propagators joining n = 4 external legs. There are two ways of joining the three 
propagators with four external lines as we see in its diagrammatic representation: 



(C(k 1 )C(k 2 )c(k 3 )C(k 4 )) t 



Using the rules in Section 0] we obtain the tree level expression for the trispectrum: 

Tf oc (k 1 ,k 2 ,k3,k 4 ) = N A N B N CB N BP [n^(k 1 )n^(k 2 )n s ^(k 13 ) + llperm.] 

N A N B N C N BBP [n^(ki)n 5 ^(k 2 )n^(k 3 ) + 3 perm.] . 



(61) 



(62) 
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5 Diagrammatic rules for non-gaussian vector and scalar field per- 
turbations 

In this section we extend the rules obtained in Section U in order to include the contributions coming from 
non-gaussian correlators of the field perturbations. We include the non-gaussianity in the form of primordial 
connected fc-point correlation functions which we will call (k — l)-spectra. With this notation k = 2 is a 
propagator. As we have studied in Section |4j in the gaussian case we have only propagators and t- vertices in the 
diagrammatic perturbative series. As an example of primordial non-gaussian contributions we have the terms 
containing the 3-point correlators of the field perturbations, the 2-spectrum (or bispectrum) B^bc m Eq. (|32|) . 
In this case, aside of the propagator and the i-vertex we have to add the bispectrum of the field perturbations to 
the building block list. The diagrammatic representation of the bispectrum in the field perturbations is shown in 
Fig. [5j The rules to evaluate perturbatively a n-point correlation function of £, including non-gaussianity in the 

(2 I :) 3 8(q l23 )B ABC (2 Jt ) 3 5(q l23 )B |JK [2nf 8(q ia ) B |Jk ° {2nf 6(q ia ) B |Jk * (2*f S(q ia ) B„- 

A 







Figure 5: The diagrammatic representation of the bispectrum in the field perturbations. Incoming momenta are 
considered to be positive. In the diagrammatic decomposition, for each scalar index I we draw a dashed line 
while for each vector index a we draw a curly line. For each vector index a we understand that there is implicitly 
an index i corresponding to the components of the vector field. 

field perturbations, are constructed by adding the (k — l)-spectra to the building block list with the appropriate 
rules to deal with their symmetries in a diagram. They are listed as follows: 

1. Draw all distinct diagrams with n-external lines, and the appropriate number s of (k — l)-spectra. The 
case k = 2 corresponds to a propagator. The order of a (k — l)-spectrum in the perturbative expansion is 
k — 1, this is (8$ £ ■ ■ ■ <5$^ ) ~ (D(P k ~ 1 ). External lines are represented by solid lines while internal lines 
are represented by dashed lines. Every vertex cannot have attached more than one external line. Every 
line in a (k — l)-spectra must be attached on at least one end to an external line, if it is attached to external 
lines in both ends it is a propagator. 

2. Label the external lines with incoming momenta k; for i = 1, ..., n and label the (k— l)-spectra with internal 
momenta pt for k = 1, ■■■,p. Label each end of each (k — l)-spectra with a field index: A,B,--- ,C. See 
Figs. [Hand El 

3. Assign the factor n^^k) to each propagator in the diagram where k is the momentum associated to the 
propagator which must flow from the end with the label A to the end with the label B (see Fig. [TJ. The 
direction of the momentum flow in the diagram is relevant when the action is not invariant under the parity 
transformation. If the momentum flows from the end with the label B to the end with the label A, the 
assigned propagator must be U§^(k) = IT^ (— k). 

4. Assign the factor Fj ^ ...^ (k!,k 2 , • • • , kfc)(27r) 3 i5(k 1 2...fc) to each primordial (k — l)-spectrum in the dia- 
gram, ki, k 2 , • • • kfc are the incoming momenta associated to the (k — l)-spectrum. See for example Fig. [5] 
for the case of the bispectrum. The Dirac delta function ensures that the momentum is conserved. 

5. Assign a factor N^ £ ...^ (27r) 3 (5(k i — pi.„t) to each vertex, a t-vertex, where the number t of derivatives 
of N is the number of lines attached to this vertex. We use the convention that incoming momentum is 
positive. The Dirac delta function ensures that the momentum is conserved. See Fig. [2] 

6. Integrate over the propagator momenta ( 2 \u J d Pi- The first n — 1 integrals can be done immediately 
using the Dirac delta functions but any further integral in general cannot be performed analytically. This 
is the case when there are integrals corresponding to loop corrections. 
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7. Divide by the appropriate numerical factor: 

• /! if there are I propagators attached to the same vertices at both ends. 

• 2 ; Z! if there are I propagators dressing a vertex. 

• p\ if there are p (k — l)-spectra attached to the same vertices at both ends, and with just one line 
attached to each vertex. 

• m\n\... if there are m lines of just one [k — l)-spectrum attached to the same vertex V\, and n lines 
of the same (k — l)-spectrum attached to the same vertex V2, and so on. 

• p\(m\n\...) p if there are m lines of each of the p (k — l)-spectra attached to the same vertex Vy, and 
n lines of each of the same p (k — l)-spectra attached to the same vertex V2, and so on. 

• l\((k — 1)\) 1 if there are I (k — l)-spectra all of the same order dressing a vertex. 

8. Add all permutations of the diagrams corresponding to all the distinct ways to relabel the k^ attached to 
the external lines. The number of permutations depends on the symmetries of the diagram: a diagram 
which is totally symmetric with respect to external lines has only one term; in contrast, a diagram without 
symmetries with respect to external lines has n\ permutations. 

In this case, given that the (k — l)-spectra are part of the building blocks of the diagrammatic representation, 
it is instructive to know the maximum number of them that we can have in every Hoop term. This information is 
relevant in order to construct all the possible diagrams at each order. In the process of drawing a Hoop diagram 
we have to take into account that the sum of the order of all the (k — l)-spectra in a diagram must coincide with 
the order of the diagram. In Table 1 we summarize the maximum number of (k — l)-spectra that could appear 
in any diagram corresponding to a /-loop correction in the perturbative series for a n-point correlation function. 



Table 1: The number of (k — l)-spectra that we need to construct the /-loop term of a n-point correlation 
function. In the rows we have the /-loop order correction and in the columns we show the maximum number of 
(k — l)-spectra that could appear at each order, [a;] represents the integer part of x. 



loop order 


spectrum 


2-spectrum 


3-spectrum 




(k — l)-spectrum 


tree: n — 1 


n-1 


[(n-l)/2] 


[(n - l)/3] 




[(n-l)/(*-l)] 


l-loop: n 


n 


[n/2] 


[n/3] 




[n/(k - 1)] 


2-loops: n + 1 


n+1 


l(n + l)/2] 


[(» + l)/3] 




[(n+l)/(fc-l)] 














/-loops: n + / — 1 


n + l—1 


[{n + l- l)/2] 


[(n + l-l)/3] 




[(n + l-l)/(k-l)] 



5.1 Examples 

Again, we are going to use the example of a single scalar field and a single vector field studied in Section [4] We 
will work out in detail the power spectrum and describe briefly the bispectrum and trispectrum. 



5.1.1 The tree level and one-loop contributions to the spectrum P^ 

For the power spectrum we have n = 2 external lines, and the tree level approximation has p = 2 — 1 = 1 
propagator, i.e. it is of order 0(P 1 ). Higher order primordial (k — l)-spectra do not appear at this order because 
the lowest higher-order correlator is the bispectrum which is of order 0(P 2 ). Hence, the tree level power spectrum 
is identical to the one in Eq. (|52|) . The l-loop power spectrum is of order 0(P 2 ) so, according to Table 1, we 
have two propagators and one bispectrum. We do not have higher-order correlators of the field perturbations at 
this order. The diagrammatic representation of the l-loop power spectrum is: 



(C(ki)C(k 2 )) 



1 — loop 






(63) 
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Using the diagrammatic rules we find the expression for the 1-loop power spectrum: 



P. 



1— loop 



c 



ki) = \n as N cd 



2 

\n a n BCDj {2n) 



J ^n AC (k 1 + p)u SB (-p) 



(2 
d 3 p 



i NANSC .l (2n) 



J ^T3n^(-kx)nc fl ( P ) + \n a n bcd J ^n^(ki)n^( P ) + 

f d 3 p 1 f d 3 p 

J 7^ B A§c( k l,P,-k 1 -p) + -N A N SC J —-B AS c(-kl,P,kl-p) 



(64) 



Comparing with the gaussian case in Eq. (|54j) . the 1-loop power spectrum has the terms in the last line as 
extra terms. These terms have a bispectrum of the field perturbations connected with a 1-vertex and a 2-vertex. 
Comparing again this expression with Eq. (|36p . we see that the one- loop contribution to the power spectrum 
of C calculated using the usual approach, employing the Wick's theorem, is fully reproduced, but this time with 
much less effort. 



5.1.2 The tree level and one-loop contributions to the bispectrum Bq 

According to Table 1, for the tree level 3-point correlation function we have a maximum of 2 propagators and 
one bispectrum of the field perturbations. With them, we can construct the following diagrams: 



(C(k 1 )C(k 2 )C(k 3 )) t 



(65) 



Comparing with the gaussian expression, we have a new term coming from the bispectrum of the field perturba- 
tions. Using the diagrammatic rules we obtain the expression for the tree level bispectrum: 

Bftkx.ka.ka) = N A N B N CD \n AC (k 1 )n BD (k 2 ) + 2pcrm]+N A N B N c B ABC (k u k 2 M)- (66) 



The one-loop contribution is constructed out of the possible diagrams that we can draw with a maximum of 
p = 3 propagators, b = 2 bispectra and t = 1 trispectrum. The resulting diagrams are shown here: 



(C(k 1 )C(k 2 )C(k 3 )> 



1 — loop 



/ 1 

• 
/ 1 


+ 


/ f 

Vr-' 

\ 
\ 
\ 








\ 1 
\ 1 




\ 
















1 / 






/ 1 

1 / 



/ 


\ 




\ 

\ 
\ 






1 
/ | 




/ 


~s 

) 




,-" 


/ . 




/ 






i I 
i ' 


\ 1 
1 




\ 







(67) 



We can see that, compared with the gaussian case, there are five new terms related to the bispectrum and trispec- 
trum of the field perturbations. Again, using the rules and doing the trivial Dirac delta function integrations, 
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the above diagrams lead to the following expressions: 



B c 1 - loop (k 1 ,k 2 ,k 3 ) = N AS N CD N BP J-^ U AC (-p)U ss (k 1+ p)U 5P (k 2 -p 



"P) + 



NaNbcNdef 



NaN S Ncdef 



NaNbcNdef 



f d 3 p 
J T^P n ^( kl ) n ss(p) n ci?( k 2 -P) + 5 perm. 



(27T) 

d 3 p 



/(^3 n^(k 1 )n S5 (k 2 )n^(p) + 2pcrm. 



(27T) 

d 3 p 



J 7^3 n^(ki)nc. £ ,(ki2)n i 3 i? (p) + 5perm. 



NaNbNcd 



f d 3 p 

I 7^T3 T n ABcn( k i> k 2,P,k 3 -p) + 2pcrm. 



-N a NbcNde 



NaNbcNqe 



\n a n b n Cde 



f d 6 p 
/ (2~^j3 n ^B( kl )- B C5£(- k 3,P,k 3 -p) + 5perm. 

/ (2^P B ^-B5(ki,P, -ki -p)n (5j g(k 2 -p) + 2pcr 

~ f d 3 p 
/ (2ttW B ^c( k i 1 k 2,k3)n 5 B(p)+2pcrm. + 

/ ?27r)3 n ^c( kl )- B s£ | s( k 2,P,-k 2 -p) + 5perm. 



-N a N b N CdS 

From the equation above, we see that the terms from fifth to nineth carry the non-gaussian contributions encoded 
in the bispectrum Babc an d the trispectrum Tabcd °f the field perturbations. We can check the numerical 
factors in front of every term by invoking the rules number 7 and number 8. 



6 Vertex renormalization 

6.1 Gaussian fields 

We can reduce significantly the number of diagrams by adopting the vertex renormalization procedure presented 
in Ref. [44]. While the discussion in Rcf. [44] is restricted to the context of multi-scalar fields perturbations, here, 
we extend the method to account for the case in which there are multi-scalar and multi- vector field perturbations. 
We describe here how this procedure works in the presence of gaussian fields and we leave the discussion of the 
non-gaussian case for next section. In this procedure the derivatives N abc-- are redefined in such a way that 
they relate the values of TV at a given background $o to the values of N at a general point x in real space. The 
derivatives of N for the background 3>o are N A g...\$ = N A §... an( t ar e related to the number of e-foldings N at 
any point x by 



N = TV($(x)) = N + N a S^>a + \n as 6Q a 6$ b + yN ABd d$ A 6$ B S$c 



(69) 



Taking the derivatives of this equation with respect to $^, then going to momentum space and taking the 
expectation value, we redefine the derivatives of N in the form 



(N 



A—Bl 



N 



1 



A-B 



-TV 



A— BCD 



d 3 p 1 

(2Tr)3~^CD(P) + Z N A-BCDEF 



d 3 pi d 3 p 2 



n CJ5(Pl) n BF(P2) 



(70) 



In the last step, we used the fact that the fields are gaussian, hence, all the expansion contains only integrals 
of propagators and not higher-order correlation functions. We see in Eq. ([70]) that the vertex renormalization 
"dresses" a vertex with an infinite series of the contractions of the vertices and the propagators integrated along 
independent internal momenta. Diagrammatically, this is equivalent to attach to the vertex an infinite series of 
loop integrals ("bubbles"). Graphically, we represent a "dressed" vertex with a solid dot as shown in Fig. [51 In 
this figure we show the expansion up to two-loops of the 3-vertex {N AB ,„ C ) . I n order to simplify our notation, 
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«• = « + ^-«? + >-< +... 

x N x . I v 

V -v V I I X 

■s N N . ' X 

X X V. X _ S N 

"V ^ -. N 

Figure 6: The diagrammatic representation of the renormalized 3-vertex N abc f° r gaussian fields. 

in the following we omit the brackets in the renormalized derivatives, then TV 'ab-G — (^AS—o)- When we 
introduce the renormalized vertex we just keep the diagrams present in the perturbative series of the correlation 
functions which do not have "bubbles" around the vertices since all the "bubbles" are already included in the 
renormalized vertex. In Fig. [7] we show an example of the type of diagrams which are allowed and the ones 
which are not allowed in the perturbative expansion of the power spectrum P<;. 

Allowed Not allowed 



+ ---*- 



Figure 7: An example of the type of diagrams allowed and not allowed in the perturbative expansion of Pf . 

6.2 Renormalized gaussian diagrams 

Taking into account the renormalization of the vertices, we can also write a set of rules for the renormalized 
diagrams as follows 

1. Draw all distinct diagrams with n-external lines and p propagators. External lines are represented by solid 
lines while internal lines are represented by dashed lines. Every renormalized vertex connect an external 
line to at least 1 propagator and is represented with a solid dot. The propagators cannot have both ends 
attached to the same renormalized vertex. See Fig. [7] 

2. Label the external lines with incoming momenta k; for i = 1, ...,n and label the propagator with internal 
momenta p^ for k = 1, ...,p. Label each end of each propagator with a field index: A, B, ■ ■ ■ ,C. See Fig. 

El 

3. Assign the factor n^^k) to each propagator in the diagram where k is the momentum associated to the 
propagator which must flow from the end with the label A to the end with the label B (see Fig. [IJ. The 
direction of the momentum flow in the diagram is relevant when the action is not invariant under the parity 
transformation. If the momentum flows from the end with the label B to the end with the label A, the 
assigned propagator must be Hs^k) = II^(— k). 

4. Assign a factor N& £ & (27r) 3 <S(ki — Pi-.-t) to each vertex, a i-vertex, where the number t of derivatives of 
N is the number of propagators attached to this vertex. We use the convention that incoming momentum 
is positive. The Dirac delta function ensures that the momentum is conserved. See Fig. [2j 

5. Integrate over the propagator momenta , 27T \3 J d 3 pi. The first n — 1 integrals can be done immediately 
using the Dirac delta functions but any further integral in general cannot be performed analytically. This 
is the case when there are integrals corresponding to loop corrections. 

6. Divide by the appropriate numerical factor: l\ if there are I propagators attached to the same vertices at 
both ends. 

7. Add all permutations of the diagrams corresponding to all the distinct ways to relabel the k^ attached to 
the external lines. The number of permutations depends on the symmetries of the diagram: a diagram 
which is totally symmetric with respect to external lines has only one term; in contrast, a diagram without 
symmetries with respect to external lines has n\ permutations. 
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As examples of these rules, we show the diagrams which represent the power spectrum up to two loops and 
the bispectrum up to one loop: 



(C(k,)C(k 2 )) 



(C(k 1 )C(k 2 )C(k 3 )) 



up to 2 — loops 



up to 1— loop 



I 

-k- 



(71) 
(72) 



They represent the expansions 



3 (up to 2 — loops) 



l - 



(kO = N A N B U AB (k l ) + -N AS N dB 



d 3 p 



nAc?(ki + p)ns 5 (- P ) 



a^ABC^DEF 



d 3 Pi d 3 p 2 



D 



(up to 1— loop) 



(ki,k 2 ,k 3 



N A N S N CB 



N AB Nc D N EP 



\n a N bc N dBP 



(2tt) 3 (2tt) 3 
n Ac( k i) n s£)( k 2) + 2 perm 
d 3 p 



n J 4D(pi) n BB(P2)nc^(ki - pi 2 ) , 



(73) 



(2tt) 3 
d 3 p 



(27T) 



^Ac(-p)^be(^i + P)lWk 2 - P) + 

i n j 4 5 (ki)nss(p)n (5i? (k2-p) + 5perm. 



(74) 



From the last equations, we see that the number of diagrams at each order reduces significantly, for instance 
the up to 2-loop power spectrum has only three diagrams instead of the seven diagrams that we have without 
renormalized vertices. In particular, the diagrammatic representation of the power spectrum is very simple 
because we have only one diagram at every order in loops. The Hoop diagram of the power spectrum is 
represented by I + 1 propagators joining two (I + l)-rcnormalizcd vertices. 

6.3 Non-gaussian fields 

We start again from Eq. (|69[) . we take the derivatives of this equation with respect to $_£, then go to momentum 
space, and finally take the expectation value. Thus, we redefine the derivatives of N in the form 



(N 



A — Bl 



= N 



A—B 



In 



A — BCD 



In 



(2tt 
d 3 pi d 3 p 2 



d 3 p 1 

— 5- n cfl(P) + T\ N A~BCDE 



3! 



n Cf>(Pl) n BF(P2) 



d 3 pi d 3 p 2 

(2tt) 3 (2tt) 3 



B, 



CDE 



(Pi 



-P2,P2-P1, 



A-BCDEF J Tj^T^p^CmPl^BFU^; "I""- . (75) 

In this case, we consider that the fields are non-gaussian so that the integrals in the expansion contain propagators 
and higher order correlation functions as well. In Eq. (|75jl we have now that the renormalized vertices will contain 
not only propagator loops but in general they will have (k — l)-spectra with all their lines wrapped around the 
vertices. Graphically, we represent a "dressed" vertex with non-gaussian fields with a solid dot as shown in Fig. 
[5] In this figure we show the expansion up to two-loops of the 3-vertex N ab ...q. We see for instance in the third 
term of such a figure that there is a 2-spectrum wrapping around the vertex. 



1 ! I 



1 ) ■ 



Figure 8: The diagrammatic representation of the renormalized 3-vertex N AB( j for non-gaussian fields. 

The diagrammatic rules for this situation are the same as the ones listed in Section [51 but excluding the 
diagrams containing (k — l)-spectra with all their lines attached to the same vertex. A general proof that the 
redefinition of the derivatives given in Eq. ([75)1 remove all the diagrams with dressed vertices was given in the 
appendix of Ref . [H] . 
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7 Conclusions 

Formulating the n-point correlators of the primordial curvature perturbation £ in terms of integrals involving 
correlators of the field perturbations is a very direct procedure; however it is in general very time-consuming, 
clutter, and not so intuitive. The same situation happened with the calculation of scattering amplitudes in 
cannonical Quantum Field Theory until Feynman exposed for the first time his diagrammatic rules (40J [41] [42] 
which, since then, have become the standard way of doing such a kind of calculations |43] ; the reason: Feynman 
rules are vey easy to remember, very easy to implement, very intuitive, and allow us to understand the physics 
behind a quantum process. A few years ago, the "Feynman-like rules" to calculate n-point correlators of £ were 
developed |44| . with the same aims as the Feynman rules in the particle physics context. In that work, only scalar 
fields were introduced as the generators of (. Since then, not many cosmologists have used this diagrammatic 
approach, wasting the interesting properties stated above. As mentioned in the introduction, the relevance of 
having an efficient method to calculate loop corrections relies mainly in the possibility that observational data 
coming from high precision cosmological probes could open a window to access effects related to such corrections; 
there are also scenarios in which evaluating loop contributions is essential because those terms dominate over the 
tree level ones and constitute a source for large and observable non-gaussianities [50] EH [EH (SHI HH ES] ■ in the 
present work, we have been interested in promoting even more this methodology among the cosmology community. 
In addition, we have extended it to the interesting case where multiple scalar and vector fields contribute to 
the generation of £; in this situation it is possible to obtain prolonged stages of anisotropic expansion (see for 
instance Ref. [53] [76] [77] [78]), as well as observable levels of statistical anisotropy [4] (20] [65] and anisotropic 
non-gaussianity [H][M][H5] which can serve as discriminators among different inflationary models PQHS]- Indeed, 
consistency relations among the different levels of non-gaussianity and the level of statistical anisotropy may be 
obtained for different classes of models so that they may easily be ruled out by observation [60] [61] [62] [79] . 
Two of the most notorious differences with respect to the multi-scalar field case of Ref. [44] are the possibility 
of parity violating interactions in the action and the possibility of anisotropic expansion; in the former case, 
the direction of the momentum flow in the Feynman-like diagrams is relevant, something which is not present 
for scalar and vector boson propagators in the Feynman rules of the particle physics context; in the latter, the 
particle production process is statistically anisotropic, which render the scalar-scalar field perturbation spectra 
and the vector- vector field perturbation spectra for the different polarizations wavevector-dependent, as well as 
making the correlators between scalar and vector field perturbations different to zero (producing scalar- vector 
field perturbation spectra) . In view of recent relevant works where the anisotropic expansion is quite prolonged 
[2"6l, |2~T1 12"51 12"91 I5B"] , we are urged to comprehensively study the generation of statistical anisotropy in that kind 
of models. 
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